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We reexamine the derivation of relativistic transport equations for fermions when 
conserving the most general spinor structure of the interaction and Green function. 
Such an extension of the formalism is needed when dealing with e.g. spin-polarized 
nuclear matter or non-parity conserving interactions. It is shown that some earlier 
derivations can lead to an incomplete description of the evolution of the system 
even in the case of parity-conserving, spin-saturated systems. The concepts of 
kinetic equation and mass shell condition have to be extended, in particular both 
of them acquire a non trivial spinor structure which describe a rich polarization 
dynamics. 



1 Introduction 

We investigate how the spinor structure of the distribution function af- 
fects the transport properties of relativistic fermions. The generahzed Wigner 
distribution function is the Green's function given by 

Fix,p)= J d^e-'P"- <^{x + ^)-9{x-^)> . (1) 

It can be developed on the basis of the Dirac gamma matrices as 

F{x,p) = fl + f^-f^ + tf^.a^'' + f^^^^^^ + i/575. (2) 

Its 16 components generally obey coupled transport equations (see e.g. 

The usual Boltzmann equation for a scalar distribution appears only if drastic 

assumptions are made. 

The spin degrees of freedom in the context of nuclear matter, for example, 
might not have been given all the attention they deserve. Relying on the 
"spin saturation" argument, the derivation of transport theories for numerical 
simulations of heavy ion collisions generally assumes that the Wigner function 
is of the simpler form F(x,p) = {j'^.Pfj, + m)ip{x,p), which permits to recover 
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the Boltzmann-Uehling-Uhlenbeck (BUU) kinetic equation. In this way, the 
only manifestation of the fermionic nature of the nucleons comes from the use 
of a Fermi-Dirac distribution function in the scalar ip[x,p). 

While the usual BUU approach is already satisfactory e.g. for the descrip- 
tion of the gross features of heavy ion collisions, and has provided valuable 
insights in the interpretation of experimental data, we nevertheless believe 
that the spinor structure should be examined in more detail. Our motivations 
are fourfold: p. 

(i) On the basis of a simplified relaxation time picturcl, a closed system 
of coupled equations for the scalar distribution / and a spin distribution /s^ 
was derived, and it was suggested that the nuclear plasma formed in a collision 
could relax more efficiently by exploring intermediate polarized states. On the 
other hand, the improvements brought lately in highly polarized beams and 
targets could be used to create a polarized excitation, and bring a new wealth 
of information on the dynamics of dense fermionic systems. 

(ii) It is interesting to remember at this stage the analogy with the case of 
the spin-polarized helium fluids, as they provide an experimental check of the 
validity of Ihe nonrelativistic limit of a transport theory with internal degrees 
of freedom^. There, effects such as increase of the viscosity, spin-heat coupling, 
spin waves in the non-degenerate regime can be observed and their magnitude 
can be compared to predictions. Similar properties can be calculated in the 
nuclear systemS. 

(Hi) The interest in such studies has been revived in the context of the 
study of thepChiral transition, using the Nambu-Jona-Lasinio as a toy model 
for the QCDEjq. In the presence of a pseudoscalar condensate, the assumption 
F{x,p) = {j'^.Pfi + m)ip(x,p) cannot be made, instead a coupled transport 
equation for the system /, /s^ is obtained. 

(iv) Besides the transport equation, the relativistic formalism provides 
also a constraint or mass shell equation. In the simple cases where the sfjiri 
saturation hypothesis is made, and also in the chiral model considered byEJu, 
this equation takes the trivial form p^ — m} . In the general case however, the 
constraint equation also retains a spinor structure. One should therefore be 
careful to use this constraint and not impose from outside the usual p^ ~ m? 
condition for all components of the Wigner function. 

The derivation of the transport equations can proceed with the standard 
Green function formalism along the Keldysh contour. In a first step, we will 
assume that the particles are only interacting via the self-consistent mean-field 
T^MpiX) — S. The transport equations for the Wigner function reduce to |^: 



"At this level of approximations, we have the properties: 'y^F^'yo = F and 7*^2^70 = S 
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[(V2)7-9 + ij-p - m) - ^e^]Fix,p) ^ 
F{x,p)[{i/2)j.d- (7.p- m) +6"^!]] 



(3) 
(4) 



Here, A is the differential operator A = d x - d p. This provides a complete 
system of differential equations describing the evolution of p) . In principle 
we could use equations (3,4) as they are. We would like however to bring them 
in a form which we could readily interpret physically as the Vlasov kinetic 
equation plus a mass shell constraint. 

The mass shell constraint is often considered as obvious, to the extent 
that some authors simply introduce it by hand at the end of the derivation 
of the kinetic equation. We would like to point out in this contribution that 
the mass shell equation is in general not a trivial one, but contains important 
dynamical and spinor information. In a first part, we show on the example of 
a simple model with scalar and vector mean fields, how different approaches 
could lead to different kinetic equations if the mass shell constraint is not taken 
consistently into account. In a second part, we compare the solution of the 
stationary problem for pseudoscalar or pseudovector mean fields respectively, 
and illustrate on these examples that the mass shell constraint acquires a non 
trivial spinor structure. 

2 The (T-u! model of nuclear interactions and the mass term issue 

We first consider the QHD-I model of the nuclear interactioiB where we have 
scalar and vector mean fields only E = (S^ + I]„p7''). 

2.1 Vlasov equation from the "quadratic" method 

A methocLto extract "kinetic" and "constraint" equations is that used by e.g. 
Elze et alQ: The equations (3,4) are put in the form 

{j.K - M)F ^0 ; F{-f.K'< - M^) = (5) 

with M = M + iN, M ^ E^, N = i9„(E,)V", if'^ = + ^V, 

= + E^;, = - 5„E,^V", V" = d/dp^. By multiplying on the 
left hand side by (7. A' + Ai) and on the right hand side by {j.K^ + A4^), and 
taking the sum and difference of the resulting equations, one obtains in the 
first order of the gradient expansion |^, after some algebra 

''Terms such as V.V/4 can be neglected. The same argument can be given performing a 
semiclassical expansion in powers of h 
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where JF^^ = 9^(1]^) — 9^(S(J). In this way, famiUar looking mass sheU condi- 
tion and kinetic-equation are recovered. This "quadratic" method however has 
been criticized lj on the ground that information could be lost in the process, 
in the same way as information is lost on the fcrmionic nature of the particle 
by going from the Dirac to the Klein Gordon equation. The issue is important 
since these authors argue that an other - "non-quadratic" - method of deriva- 
tion yields a different kinetic equation where the mass term Af9Q,(Ss)V" is 
absent. We will return to this issue in the next section. We first analyze more 
carefully the equations and show that the mass term is also obtained without 
squaring the Dirac equation. 

2.2 A "pedestrian" method 

The components of system of equations (3,4) on the basis of Dirac gamma 
matrices yield in the first order gradient approximation 

V^r + dx{M)W^f^O (8) 
T^J + 4n"/„^ + dx{M)V^U = (9) 

_ igM.«/3p^ + n^'^r^ - dx{M)^^r = (10) 

^^"''^V^f^p - 2Wh + dx{M)V^J^ = (11) 

2n^f^ + dx{M)^^h = Q (12) 

li^r - Mf = (13) 

n^/ - M/^ + VaF" ^ (14) 

T^l^M - ^Mf^, + 2e^,„^n"/5'3 = (15) 

2?m/5 - 2M/5^ + 2e^,^0l]!'r^ = (16) 

V^f^ + 2M/5 = (17) 

where a'^fc'^l = a^^b^ — a^h^^, et^'^"^ is the completely antisymmetrized Levi- 
Civita tensor and we have defined as before 11^ = + S^, AI = m — J^s, 
= d/dpx, V^ = d^~ ^aH^V^ 

In order to calculate the physical observables (baryon number current 
and energy momentum T'^'^), only the components / = Tr{F) and = 
Tr{'j'-'-F) are necessary. The component allows to calculate the nucleon 
and vector field contributions to and T^" while / gives the gap equation 
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and the scalar field contributions We can obtain a kinetic- like equation by 
multiplying (9) with 11^: 

n^a^/ - d^{M)r - [u^dxin^) - Mdx(M)] vV - o. (is) 

A kinetic equation for can be derived from (9) and (15) as follows. We 
multiply (9) with M and (15) by H'^ and sum the resulting two equations in 
order to eliminate the f'^'^ contribution. We obtain in this way 

n„(9"/'3-a'^(M)/-[n„aA(n")-M9A(M)]vV'^+[a''n"-a"n'3]/„ =0.(19) 

Equations (18) and (19) can be rewritten as 

[n"2?„ - Af9„(E,)V"] / + a,,(S,)/^ = (20) 

[n"Vo. - M9„(S,)V"] r + ^^Va = (21) 

It is straightforward to see that these equations are exactly what one obtains 
from taking the trace of (6) with I and 7'^ respectively. As for the mass shell 
equations, we obtain the one for / by multiplying (14) by 11'' and using (13) 

(tf - M^)f + n^P,/^" = (22) 
and the one for by multiplying (10) with If making use of (13) and (14) 

{n^~M^)r +(i/2)e''"'3^n„i?^/5A + Afi?„/''" - aA(s,)n,vV^'^=o. (23) 

We notice that (22,23) differs from the traces of (7) with / and j'^. 

We now proceed to show how can one recover the equations of Elze. Let 

us apply derivatives^ to (9), obtaining n^V^f" = df,{Y.y^)f'"' + and 

dxi^s)TluV^F'' = -dxii:,)F^ + 0^; and to (15), obtaining e'^^^^H^P^/sA + 

2MVaf^"^ = e''"^^X>„(n;3)/5A + O^. Substituting in (22,23) we see that the 

traces with I and 7'' of the mass shell equation (7) of Elze et al. are indeed 

recovered. The same procedure could be applied to obtain kinetic and mass 

shell equations for all components of the Wigner function. 

Note that (20-23), can also be obtained directly by multiplying (3,4) right and 

left by (7.n + M)0 and take the sum and difference 

2(tf-M2)i^ 

+ ^n^P4a'''', F} + '-MV^ [F, F] - ^n^9„(S,)V"[7^, F] = (24) 
[n''P^-M9„(E,)]F 

+ in^I?,K'', F] + ]^MV^ {F, F} - ln^(}„(E,)V"{7'', F} = Q. (25) 

■^Only at this stage an operation is made that cannot be inverted 

''without derivative, one can go back any time to (3,4) multiplying again with {7.11 — M) 
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Then, (22) and (23) arc the traces of (24) with / and 7^. 
Until (24,25), we can be fuUy sure that no information was lost. It is tempting 
to use (6,7) since these equations take a more compact form and have the nice 
feature of containing explicitly the precession terms J^^i, [cr^" , F] . They will 
probably give no trouble for all practical purposes. One should however keep 
in mind that the system (24,25) is more general than (6,7) as some loss of 
generality is taking place in the derivation process. In any case, we find that 
the mass termJljf9a(I]s)V" is present in (20,21,25) contrary to the claims of 



the authors of We will examine their equations in section ( p.3[ ) and show 
that the mass term discrepancy stems from a faulty mass shell condition. 

We now discuss a much more severe restriction commonly imposed on the 
mass shell, which consists of approximating (7) by (11^ — M'^)F = 0. If / and 
/'^ are remaining on the mass shell 11^ — = 0, equations (20,21) would 
represent a self contained system for the evolution of / and f^. From (22,23) 
we see that in general however they deviate from this mass shell by terms 
involving the derivatives of /^'^ and f^. Only in symmetric nuclear matter 
in a state of (homogeneous) equilibrium, and when the system is unpolarized, 
does the stationary solution reduce F = (7.n+M)(5(n^— M^)^?. Unfortunately, 
very often it is assumed that the spinor structure remains unchanged out of 
equilibrium, (F = / + /^7^, = /^"^ = /g = /g^). Then the system (8-17) 
simplifies, namely: Eqs. (11,12,16,17) contain no information. Due to the 
assumption f^"' = 0, Eq. (14) gives = (H^/A/)/ while Eq. (13) returns 
the mass shell condition (H^ - M^)/ = 0. Eq. (10) gives nl'^/'^l = which is 
trivial with ft" = (H^/M)/. We are left with Eqs. (8,9,14), and also (18,19) 
which were derived from the former. Then, either 

* by replacing f^' by (H/'/M)/ in Eq. (8), 

* or multiplying (9) by 11'^ and replacing 

* or replacing in Eq. (19) which was obtained from (9) and (15), 
we arrive at 

{n^9^ - [Wdx{lV^) - A/aA(Af)V^] }(^) = O. (26) 

This equation might be further transformed by taking the change of vari- 
ables (a;,p) {x^Ti) so that we obtain 

[wd, - n'^^^^^ - Md,{M)^]}{l^) = 0. (27) 

with = d^^{Yj^) — d^{Y.'^). This is the equation which is extensively used 
in numerical simulations of heavy ion collisions with a Boltzmann-Uehling- 
Uhlenbeck collision term on the right hand side. 
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Since Eq. ( p7| ) was arrived at by three different methods, we could think that 
the result is indeed safe. However we still not have used all the information 
available in the system (8-17). For example, we could have replaced by 
ni^/Mf and take f^"' = in Eq. (15) directly and obtain 

"'"^"'(^^ ~ n[^9A(n''i)v^(^) = 0. (28) 

Strictly speaking, we simply do not have a justification to put — 0. This 
term has no reason to vanish out of equilibrium. It could for example ac- 
quire contributions like J-^'^f, IlI^S'^l (A/)/, or time-reversal invariance break- 
ing terms like U^f^f^. 

2.3 Vlasov equation from the "Zubarev" method 

In this alternative derivation ("Zubarev formalism"), after defining a 
proper time r = x^u^^ one puts the Dirac equation into Hamiltonian form 

dip -iA 
iu'^dn.ip = = Hip with H = —ia^^Undjj -\- m^.u + ^.uT,e ^ 

OT 

The conjugate equation for ip reads 

iu'^df,'^ = ^^H withH = i(7'"'uf,K~m'y.u-e'^Y:j.u. (29) 

OT 

We have H = — 7oi?^7o and E = 70X1^70 = S. Acting on the definition of F 
with u.dF — ... and making these replacements, we arrive at^ 

u.dF ^ J d'^Re-'P '^ -ilpH ® iP - ii> ® Hijj = -i{HF + FH). (30) 



Choosing pDw m'' = II^/v 11^, one arrives at the "kinetic equation" of Smolyan- 
sky et al. H. 

li.dF + in^9,K^ F] + jM[7-n, F] - n^a„(5],^)v«i^ 

-in^9„(I]„,) K^ V"F] - \d^{Y.s) {7.n, V"F} = O (31) 

One can easily see that the same equation can be obtained by multiplying the 
eqs. (3,4) by 7. 11 from the left and from the right and taking the sum. 

However, we still have an other equation, obtained by multiplying on the 
right and left hand side and subtracting, which is the corresponding constraint 
equation. It reads 

'^Note that this can be put in the form of the Von Neumann equation idoF = [H, F] if we 
choose the coordinate along u*^, so that = and H = —H in the mean field 

approximation. 



7 



§n^9^ {cr'^'', F} + 2U^F - M {7.n, F} 

-§^^a<,(I],,){a''^v"^^}- i[7.n,a„(E,)v"(F)] = (32) 

Until now, the system ( ^Ip^ ) is fully equivalent to the original set of equations 
(3,4)0. However, the second "constraint" equation is not taken into account 
by Smolyansky et al. who replace it by the mass shell condition li^ = M^, 
although we can see from the explicit expression that the A/^ term is missing 
and that it contains instead dynamics at the same order as the other equation. 
Moreover, one can see that in the free case, the well known kinetic equation 
p.dF = 0, is not recovered from (|3^) alone ^ This example illustrates that 
some care has to be paid that the seemingly obvious mass shell condition 
should be obtained from the calculation, less some part of the dynamics may 
be neglected by imposing this condition beforehand. 

2.4 Vlasov equation: discussion 

There are several ways one can derive an equation of the form Il.dF+DK[F] = 
0, each method of derivation leading to a different set of additional terms 
Dk[F]. To this equation should be associated the corresponding mass shell 
condition which also contains dynamical effects {iP — M'^)F + Dm[F] = 0. If 
the usual mass shell condition were imposed in the form (H^ — M'^)F = from 
outside, terms Dm[F] of the same order as Dk[F] would be neglected, leading 
to a wrong description of the dynamics. The Zubarev method is particularly 
treacherous in this respect, since the neglected terms are of order zero in the 
gradient expansion. The method developed by Elze et al. on the other hand 
permits to minimize the importance of the Dm[F] terms so that these are 
of second order in the gradient expansion if the interactions are switched off. 
When the mean field is present, however, Dm[F] remains of first order and has 
to be considered. 

3 Pseudoscalar vs. Pseudovector mean fields 

In order to illustrate the highly nontrivial structure of the mass shell 
constraint in the general case, we will compare here the stationary, homo- 
geneous solution in the presence of a nonvanishing pseudoscalar or pseudovec- 

f Since we can always recover the initial set of equations by multiplying a second time with 
7.n and use (■y.Tl)^ = 

9It can be shown that the equations p.dF = and (p^ — m?)F = can be recovered 
from the system (31), (32) by submitting it to a series of manipulations which amounts to 
return back to the system (3,4) by multiplying (31) and (32) by 7-11, and taking the sum 
and difference, and then apply the technique exposed in the previous section 
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tor mean field. The equations (3,4) now read (7.73 — m — T,)F{x,'p) = = 
F{x,p){'^.p — TO — E) The sohition of these equations is found to bcQ 

Fix,p) = Pl{x,p) Q{x,p) Pr{x,p) 5{D) (33) 

Here, Pl{x,p) and Pu{x,p) are matrices in the Chfford algebra such that 
{■j.p — TO — S) Pii{x,p) = D = Pl(x,p) {-j-p — to — E), D = det{j.p — to — E) 
being a c-number. The Dirac S function provides for the cancelation of the 
products (7.P — m — E)_F and F{'y.p — m — E). Q{x,p) is an arbitrary matrix 
in the Clifford algebra and contains the thermodynamic averages. It can be 
shown that Pr = Pl = P, and P can be interpreted as a projector on the 
positive-energy state on the mass shell D — 0. 

3.1 Scalar, vector and pseudoscalar mean fields 

We assume here that the mean field self energy is of the form: E = (Eg + 
'^vfj.j'^ + S575). The stationary homogeneous solution to the Vlasov system of 
equations is now found to berj 



{pl-ml-j:l)f^O ; F = ^f ; h^^f 

TO* ?71» 



with again p'^ = p^ + E"^, to* = m — E^. The solution factorizes as 

F{x,p) = Pe [A©<5(p^ -ml- ^iMx.p)] P^ (35) 

with generalized definitions of two projection operators: The projector on the 
positive energy state of the mass shell pi — ml — Eg = 

^^e = 1 + — 7m + *— 75 (36) 

TO* TO* 

and the projector on the local spin direction 

Ae-(l+757M5^)/2; = (37) 

(m^ + E^)^ 

with / — 5{pI — ml — E|)iy9 and the constraint p^^S^ = 

In this case, we have four independent distribution functions: / and the three 

components of fl^ orthogonal to p* . 

''The reader may easily convince himself that inserting the expression for f^^'^ in that of 
leads to an identity 
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3.2 Scalar, vector and pseudovector mean fields 

Let's assume that S = (S5+S^^7''+S5^757'*). We find that the stationary 

homogeneous solution of the Vlasov equation contains only one independent 
distribution function (the scalar one /). The other components are given by 

f^^ = A^+B^ ; f^=C^+V^ ; /5 = ; /z'- = e'*-P^^^f (38) 

m* m* TO* TO, mj 

The expressions of A,B,C,T>,£ are given in the appendix. The mass shell 
condition is given by —mt,f + p^f" — Ss^/" = or equivalently Df = with 

£>=[(!- pl/ml - EsaS^/™?)' - 4(p;E^/to2)2 + A^.^J:'^ /ml] (39) 

= mt- 2ml {P*a + Ssa) (Pa " ^Sa) + (Pa " ^Sa) (?>^ " S5/3) (j^a + ^Sa) (p^ + ^5/?) 

The solution can again be expressed in a compact form with the help of a 
generalized projector Pq 

F{x,p) = Pe5{D),p{x,p) (40) 

p _ jm. -p;7^ - ^5^757^) (m. -ptY - ^5.757") (m. -p;^'' - SspTslO 
® 2m,(m2-p2 + S5aS^) 

on the mass shell D = defined as above (eq. 39). It is obvious that the 
physics described by the solution in a pseudoscalar field will be very different 
from the physics in a pseudovector field. In the case of the Walecka model for 
example, it will thus not be a simple matter of convenience to choose PS or 
PV coupling of the pion field. 

4 Conclusion 

On the basis of two examples, we have wished to stress the importance 
of taking into account the mass shell condition properly. Only in very simple 
cases does this condition take the simple form (II^ — M^)F = 0, in general 
however it is highly non trivial. 

In the case of the a - w model of nuclear matter, it was shown that one cannot 
separate clearly the "kinetic" and "mass shell" conditions in general. The mass 
shell condition always contains some dynamical contribution in the interacting 
case. Overlooking this fact may even lead to wrong dynamics. 
In the case of a system with pseudoscalar or pseudovector coupling, we have 
solved the stationary case and shown that the mass shell condition acquires a 
rich spinor structure. The PS or PV coupling lead to very different solutions. 
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Appendix 

The full expression of the coefficients entering the PV result is: 

A = [b^E^K (/ +K/7m*) - {pl/ml - l)E5„p7m^/] /T 

B = [{pl/ml - l)(S5„/57m. - /) - (S5aP«/m2)V] 

C = [(1 + S5„S«/m^) {f+p*J-/m.) - (S5„p:/m^)V] /J" 

V = [(E5oP:/m2)5]5„/5«/m* - 2/ - SgaS^/m^/)] /jr 
f = -(y^ + I?)/2 

= [(1 + E5aS;^/m2)(l + (I]5aP?/m2)2 ^ pl/ml] /{J^G) 

T = {J:5c.p"/mlf + (1 - pl/ml)il + SsoS^m^) 

= [-2(pl^^/mif+pl/ml (pl/m'i-l + i:,^J:^/ml)]/gmJ 
Q = pl/ml - 1 - SsaS^m^ 
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